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Abstract 

We discuss the effect of horizontal fluxes on the accuracy of a conventional plane-parallel radiative transfer 
calculation for a single pixel, known as the Independent Pixel Approximation (IPA) at absorbing wavelengths. 
Vertically integrated horizontal fluxes can be represented as a sum of three components; each component is the IPA 
accuracy on a pixel-by-pixel basis for reflectance, transmittance and absorptance, respectively. We show that IPA 
accuracy for reflectance always improves with more absorption, while the IPA accuracy for transmittance is less 
sensitive to the changes in absorption: with respect to the non-absorbing case, it may first deteriorate for weak 
absorption and then improve again for strongly absorbing wavelengths. IPA accuracy for absorptance always 
deteriorates with more absorption. As a result, vertically integrated horizontal fluxes, as a sum of IPA accuracies for 
reflectance, transmittance and absorptance, increase with more absorption. Finally, the question of correlations 
between horizontal fluxes, IPA uncertainties and radiative smoothing is addressed using wavenumber spectra of 
radiation fields reflected from or transmitted through fractal clouds. 
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1. Introduction 

In order to correctly interpret radiation reflected from or transmitted through clouds and 
measured by satellites and ground-based radiometers or two aircraft flying above and below 
clouds, we need to better understand interactions between inhomogeneous clouds and solar 
radiation. The discrepancies between observed and model-predicted shortwave absorption (e.g.. 
Cess et al., 1995, Wiscombe, 1995), between cloud optical depths estimated from satellites and 
ground measurements (Min and Harrison, 1996), between single scattering albedo retrieved from 
in situ radiation measurements and computed from measured droplet size distribution (Pincus et al., 
1997), amongst other reasons, are strongly affected by cloud horizontal inhomogeneity. 

Net horizontal photon transport (i. e., horizontal fluxes) are a direct consequence of the 
inhomogeneity in cloud structure. Horizontal fluxes and their effect on the accuracy of the pixel- 
by-pixel one-dimensional (ID) radiative transfer calculations has recently undergone close scrutiny 
for conservative scattering (Marshak et al., 1995a; Barker, 1996; Davis et al., 1997a, b, c; 
Chambers et al., 1997; Titov and Kasjanov, 1997; Titov, 1998; Zuidema and Evans, 1998). 
However, the properties and magnitude of horizontal fluxes in absorbing wavelengths are still 
poorly understood. As far as we are aware, only Ackerman and Cox (1981) and Titov (1998) 
briefly discussed correlations between horizontal fluxes at different wavelengths. 

This paper partly fills this gap. We discuss here an intriguing question of whether the 
accuracy of the Independent Pixel Approximation (IP A), as a ID radiative transfer approximation 
for each pixel, is better for transparent or absorbing wavelengths. The main original points to be 
addressed here are: 

• dependence of net horizontal fluxes on single-scattering albedo; 

• connection between pixel-by-pixel accuracy of the IPA and horizontal fluxes; 

• radiative smoothing and horizontal fluxes in wavelengths with liquid water absorption. 
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2. Cloud Models and Numerical Radiative Transfer Tools 

The key point of our cloud models is a scale-invariant horizontal distribution of cloud 
optical depth x. We started with simple one- or two-dimensional (2D) “bounded cascade” models 
(Cahalan, 1994, Marshak et al., 1994) that simulate cloud inner structure with x > 0 for the whole 
domain. Then we supplemented this structure with gaps using a simple procedure that first linearly 
transforms x and then sets negative values to zero (Marshak et al., 1998). Finally, we added 
cloud-top variability using a fractional Brownian motion (Mandelbrot, 1977). 

We show results where we set mean optical depth <x> = 13, mean geometrical thickness ( h ) 
= 300 m, bounded cascade model with variance parameter p = 0.3 and scaling parameter H = 1/3 
(Marshak et al., 1994). All these parameters are typical of marine stratocumulus (Cahalan et al., 
1994). Pixel size was set to 25 m; a 10-step cascade model then yields outer scale L = 25.6 km. 
All results presented below are averaged over ten independent realizations of stochastic cloud 
models. 

Both Henyey-Greenstein with asymmetry parameter g = 0.85 and Cl cloud (Deirmendjian, 
1969) scattering phase functions were used. To study the dependence on single-scattering albedo 
GJo, seven albedos were chosen for our simulations: 1.00, 0.999, 0.996, 0.99, 0.98, 0.95, 0.9 
which cover the range currently calculated for pure liquid water clouds in the shortwave spectrum. 

Two traditional numerical methods were chosen for radiative transfer calculations in the 
cloud models above. They are the Independent Pixel Approximation (IPA) and Monte Carlo (MC) 
simulation. While the later is the “exact statistical” solution (Marchuk et al., 1980) of the three- 
dimensional (3D) radiative transfer equation (see below) with periodic boundary conditions, the 
former is a plane-parallel solution applied on a pixel-by-pixel basis (Cahalan et al., 1994) and 
calculated using either a two-stream approximation (Lenoble, 1985), DISORT (Stamnes et al., 
1988) or MC with very large horizontal pixel sizes. 
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3. Horizontal Flux and Its Components 

To determine photon horizontal transport, we start with the radiative transfer equation (e.g., 
Lenoble, 1985) 


Q.WI = - c(x)7(x, Q) + a(x)G5 0 J p(ft’->Q)7(x,Q’)dQ’ ( 1 ) 

471 

where 7(x,Q) is a radiance at a point x = (x,y,z) in direction Q = (Q^Q.^), p(Q’->Q) is the 
(normalized) scattering phase function, a(x) is extinction coefficient and, finally, G3 0 is a single- 
scattering albedo. If Eq. (1) is integrated term-by-term with respect to Q. (over 4jt) and z (from 
cloud base, zb, to cloud top, z t ), we get a relationship (Titov, 1998) between vertically integrated 
horizontal fluxes H(x,y), column absorption A(x,y), reflectance R(x,y) and transmittance T(x,y). 
In a simple case of a(x) = ct(x), we have (Ackerman and Cox, 1981; Davis et al., 1997a; Marshak 
et al., 1998) 


f 5 f Zl 9 Zt 

H ( x > J 5x J &xI(x,z-,Cl)dQdz=- \ g- J Q z /(x,z;Q)dQdz-(l-GJ 0 )<J(x) f f 7(x,z;Q)dQdz 

Zb Zb 4it Zb 4k 


The first term in the right-hand side of the above equation is the difference between two net fluxes 
at z t and Zb> while the second one is column absorption. In other words, 

H{x) = { [ l-/?(x)] - [T(x)~ 0] } - A(x), 0 < x < L, (2) 

i.e., vertically integrated horizontal fluxes H(x) is determined as a difference between “true” 
column absorption A(x) and its “apparent” counterpart, 1 -/?(.«:) ~T(x). 

Equation (2), however, does not say anything about computational accuracy if, instead of 
full 3D radiative transfer, a fast IPA is used. We know that IPA treats each pixel as an independent 
plane-parallel medium neglecting any net horizontal photon transport. In order to relate horizontal 
fluxes and pixel-by-pixel accuracy of the IPA, we substi ute unity in Eq. (2) by the sum of 
#IPA(*). TipaW and Aipa(x). As a result, 77 is represented as a sum of three components, 

H(x) = Hr(x) + Hjix) + Ha(x). (3) 

Each component in Eq. (3) is a pixel-by-pixel EPA error (or accuracy). 
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Hf{x) = FipaW - F(x) (4) 

where F is either R, T or A. Since Hr, Hj and Ha are reflection, transmission, and absorption 
components of vertically integrated horizontal fluxes, we also will call them horizontal fluxes for 
photons reflected from cloud top, transmitted to cloud base, or absorbed by cloud column, 
respectively. 

To measure the magnitude of horizontal fluxes, we will use the norm, 

L 

IIFll = [ J IF(;c)l 2 dx] 1/2 (5) 

0 

where L is the outer scale or the size of the basic cloud cell. It follows from Eq. (3) and the norm 
definition (5) that 

llflll < HHrH + II// 7 H+ llff A ll. ( 6 ) 

Next we discuss the dependence of each component in Eq. (3) on single-scattering albedo, 
GJo> and examine whether or not the magnitude of horizontal fluxes is related to radiative smoothing 
(Marshak et al., 1995a) and the accuracy of IPA. We will also examine how close the left and right 
parts of ( 6 ) are for two different solar angles. 

4. Dependence on Single-Scattering Albedo 
4.1 IPA accuracy on a per-pixel basis 

Let us “coarse-grain” the above horizontal fluxes (H, Hr, Hj and Ha ) to scale r, 

x +r 

F{r r x) = | J F(x’) dx’ (0 < x < L-r, 0 < r < L). (7) 

* 

It is natural to assume that if scale r in Eq. (7) grows, the assumptions of IPA become more 
justified, the IPA itself more accurate, and the horizontal fluxes ll//f<r)ll smaller. Figures la, lb, 
and lc show that this is true for Hr, Hj and Ha, i.e., 

\\Hf(r)\\ — » 0, r — » L; F = R,T, and A. ( 8 ) 
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It follows from inequality (6) that (8) is valid for H as well. The effect of single-scattering albedo 
GJo on horizontal fluxes is, however, different for reflected, transmitted and absorbed photons. 

Figure la illustrates the dependence of Hr on both r and GJo. We see that the more 
absorption the shorter photon horizontal transport for reflected photons. As a result, the IPA pixel- 
by-pixel error decreases with the decrease of GJ 0 . Note that, for GJo=0.9, IPA is almost accurate 
even on a per-pixel basis (Marshak et al., 1995b). This is expected since the contribution of 
multiple scattering to the albedo field decreases and very few photons travel between pixels. 

The situation with transmitted photons is surprisingly different: the accuracy of the IPA on 
a per-pixel basis first decreases (from GJ 0 = 1.0 down to GJ 0 = 0.98) and only for very strongly 
absorbing wavelengths (GJo > 0.98) does it increase again (Fig. lb) reaching the accuracy level of 
the conservative scattering in case of GJ 0 = 0.9. To explain this, note that in contrast to standard 
deviation of /?ipa(x) which monotonically decreases with the decrease of GJ 0 the standard deviation 
of 7ip a (t) first increases and then decreases again for strongly absorbing wavelengths. [This can 
be shown with a little algebra for all solar zenith angles within the 2-stream approximation (e.g., 
Lenoble, 1985).] A good linear relationship (not shown here) between standard deviation of 
/?IPA(x) and Hr from one side, and standard deviation of Tip A (x) and Hj, from the other side, 
completes the explanation of both Figs, la and lb. 

The increase of pixel-by-pixel IPA absorption errors, IIHa(GJ 0 )II, with more absorption 
(Fig. lc) is understandable; it follows directly from both the natural increase of Ajpa(x) itself and 
its standard deviation with increasing of co-albedo 1— GJ 0 - 

Figures 2a and 2b illustrate the joint effect of all horizontal fluxes [see Eq. (3)] for both 
high (0o = 0°) and low (0o = 60°) Sun. The general tendencies of horizontal fluxes are similar for 
both solar angles (with the increase of co-albedo, ll// /4 (GJ 0 )ll sharply increases, ll// /? (GJ 0 )ll 
monotonically decreases, and II//t(GJo)II slowly decreases for 0o = 60° and first increases and then 
decreases for 0 O = 0°). However, the vertically integrated ho izontal fluxes II//1I is much closer to 
the sum ll//All+ll//rll+ll/// ? ll in case of slant illumination (Fig. 2b) than in case of zenith Sun (Fig. 
2a). 
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To interpret this, note that for high Sun horizontal fluxes for reflected and transmitted 
photons are mostly anticorrelated, while for low Sun they are mostly correlated. This is a direct 
consequence of radiative channeling (Davis, 1992; Davis et al., 1997, 1998). As a result, for the 
majority of pixels the product HrHj is negative if 0 O = 0° and positive if 0 O = 60°. Thus 

IWII « lltf/fll + \\Ht\\+ ll/f^ll, 6o = 0°; (9a) 

lltfll « llf//?ll + WHt\\+ llff A ||, 0 O = 60o, ( 9b) 

as we see in Figs. 2a and 2b. 

Finally and most importantly, the increase in vertically integrated horizontal fluxes ll//(GJ 0 )ll 
with the increase of co-albedo 1-GJo is entirely determined by the increase of ll//^((Do)ll. In the case 
of slant illumination, this is true for all single-scattering albedos while in case of high Sun only for 
strongly absorbing wavelengths. To conclude, horizontal flux H{x) defined in Eq. (2) is not 
directly related to IPA accuracy for either reflectance or transmittance; for strongly absorbing 
wavelengths, the IPA errors Hr{x) and Hjix) can be sufficiently small but nevertheless horizontal 
fluxes H{x) are large because of the absorptance error H\{x). 

An original explanation of the increase of vertically integrated horizontal fluxes with co- 
albedo was advanced by Titov (1998). He stated that, only in case of conservative scattering, 
photons traveling between neighboring pixels do not contribute to the increase of horizontal fluxes; 
in absorbing cases, a photon traveling back and forth horizontally changes its “weight” after each 
order of scattering. As a result, the more absorption, the larger changes in photon’s weight, the 
bigger its contribution to horizontal fluxes. 

4.2 Radiative smoothing and horizontal fluxes 

In the case of conservative scattering, Marshak et al. (1995a) defined radiative smoothing 
as a radiative transfer process that diminishes the small-scale fluctuations of cloud structure 
producing much smoother (sometimes almost everywhere differentiable) radiation fields. Radiative 
smoothing is a consequence of horizontal fluxes that are driven by the gradient in cloud structure . 
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The scale q that characterizes this process is called “radiative smoothing scale.” Based on the 
diffusion theory, Marshak et al. (1995a) estimated T| as 


T| = 


<A>[(1-S)<t>]- 1/2 , 

(h). 


for reflected radiation 
for transmitted radiation 


( 10 ) 


Here ( h ) and (x) are the cloud’s mean geometrical and mean optical thicknesses, respectively and g 
is the scattering phase function asymmetry factor. Davis et al. (1997b) show that the estimation 
(10) agrees (at least for reflected radiation) with the spatial correlations observed in Landsat scenes 
of marine stratocumulus by Cahalan and Snider (1989). 

As was noticed earlier (Marshak et al., 1995a; Davis et al., 1997a, b), the radiative 
smoothing scale T| serves as a critical value where IPA effectively breaks down: for scales smaller 
than T|, real radiation fields are much smoother than their IPA counterparts which have a similar 
variability to cloud structure. Since the radiative smoothing scale also “measures” the scale where 
horizontal fluxes become important, there are two important questions to ask: How well does T) 
characterize the magnitude of vertically integrated horizontal fluxes H and how is it related to the 
pixel-by-pixel IPA errors, Hf defined in Eqs. (3) and (4)? 

To answer these questions we consider the effect of different single-scattering albedos G5 0 
on the radiative smoothing scale q. Figures 3a and 3b illustrate wavenumber spectra E(k) of the 
nadir radiance fields for three different G5 0 = 0.999, 0.98, 0.95 and two solar zenith angles 9 0 = 0° 

(Fig. 3a) and 0 O = 60° (Fig. 3b). Both IPA and MC spectra aie plotted. We start with the simplest 
case of Sun in zenith. 

First, we see that all three IPA reflection fields have scale-invariant wavenumber spectra, 

E(k)~l/$ (ii) 

with p = 1.5, similar to the cloud optical depth field (see Fig. 3c). In contrast, all three MC fields 
exhibit a scale break that is determined by the radiative smoothing scale q and is a direct indicator 
of photon horizontal transport. Next, the scale break moves towards smaller scales with the 
decrease of G5 q> i*e., 
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Tl(GJo) > Tl(Oo’), ifGj 0 >GJ 0 \ 


( 12 ) 


where T|(GJo=l) is estimated by Eq. (10). In other words, inequality (12) states that the more 
absorption the smaller T|; this is understandable if the rescaled optical depth in case for 
conservative scattering (l-g)(x) [see Eq. (10)] is substituted with (1-GJogXt) for absorbing 
wavelength (van de Hulst, 1980, Platnick, 1997). 

Finally, we note in Fig. 3a that the smaller r| the flatter the small-scale slope; this means 
less radiative smoothing, thus more accuracy in the IPA reflection. This is in good agreement with 
the conclusion of the last section and Fig. la: the pixel-by-pixel IPA error in the reflected radiation 
monotonically decreases with GJo- From the other hand, the vertically integrated horizontal fluxes 
H defined in Eq. (2) increase while the horizontal fluxes for reflected photons Hr defined in Eqs. 
(3)-(4) decrease. Therefore, horizontal fluxes H are determined for strongly absorbing 
wavelengths mostly by the IPA absorption error Ha and are not directly related to the radiative 
smoothing scale r\ which relates to Hr, at least for zenith Sun. 

The situation with oblique illumination (Fig. 3b) is more complex. As a direct consequence 
of optical shadowing, the MC reflection fields have much more fluctuations at intermediate scales 
than their IPA counterparts (Marshak et al., 1997). The small-scale behavior is still governed by 
radiative smoothing, but this is not clearly seen in the wavenumber spectra. In summary, there are 
two competing processes that determine the wavenumber spectra: shadowing (or “roughening”) 
and smoothing. While the former one flattens the spectra, the later steepens it; roughening is more 
pronounced for intermediate scales while smoothing for small scales. 

Finally, Fig. 3c illustrates wavenumber spectra for the high Sun radiation transmitted 
through clouds. As we see, the inequality (12) is not valid anymore. Indeed, comparing Fig. 3c 
with Fig. 2a, we find that scale break and the radiative smoothing scale r\ is determined by the 
horizontal fluxes for transmitted photons Hj which does not show much variability with respect to 
GJo (Figs, lb and 2a). Another indication of weak dependence of r| on GJo is Eq. (10) for 
transmitted radiation: it depends only on mean cloud geometrical thickness ( h ). Hence, similar to 
reflected fields, radiation transmitted through clouds exhibits a scale-break which is characterized 
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by the radiative smoothing scale T| and depends only on the IPA accuracy for transmittance H T and 
not on vertically integrated horizontal fluxes H. 

Similar to Fig. 3c, wavenumber spectra for the radiation transmitted through clouds from 
the low Sun (not shown here) does not exhibit a strong G5 0 -dependence. Because of longer photon 
path, the small-scale behavior of transmittance for oblique illumination is smoother than in case of 
zenith Sun. In contrast to up-welling radiation (Fig. 3b), the shadowing effect is well seen only 
for gappy clouds. Note that because of shadowing for oblique illumination, it is very hard to 
estimate from wavenumber spectra the smoothing effect of horizontal fluxes and IPA accuracy. 
This statement is valid for both up- and down-welling radiation. 

5. Summary and Discussion 

Vertically integrated horizontal fluxes, defined as the difference between “true” and 
apparent (measured) absorption in Eq. (2), can be represented as the sum of three “horizontal 
fluxes” for reflected, transmitted and absorbed photons [see Eq. (3)]. These fluxes are also pixel- 
by-pixel IPA errors for reflectance, transmittance and absorptance, respectively. 

We show that, in general, the magnitude of vertically integrated horizontal fluxes H 
increase with the increase of single scattering co-albedo l-tD 0 ; this is true for all G5 0 in case of low 
Sun and for strongly absorbing wavelengths (1— G3 0 > 0.01) in case of high Sun. However, the 
increase of their magnitudes is not correlated with the accuracy of IPA for reflectance and 
transmittance. The accuracy of IPA on a pixel-by-pixel basis for radiation reflected from and 
transmitted through cloud is described by the net horizontal transport of reflected {Hr) or 
transmitted {Hf) photons, respectively; both are just components of the vertically integrated 
horizontal fluxes H. In contrast to H, the magnitude of Hr arid Hj decreases with the decrease of 
GJo (Hr for all solar angles and Hr for low Sun) or weakly sensitive to G5 0 (Hr for high Sun), 
while the increase of vertically integrated horizontal fluxes H is due to the increase in the error of 
IPA absorption on a pixel-by-pixel basis. 
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The radiative smoothing scale T|, different for reflected and transmitted radiation [see Eq. 
(10) and Figs. 3a and 3c), is a measure of the pixel-by-pixel accuracy of IPA. At least for high 
Sun and reflected radiation (Fig. 3a), t| decreases with more absorption [see inequality (12)]. This 
effect is more hidden for low Sun because of shadowing and thus roughening at intermediate scales 
(Fig. 3b). For transmitted photons, r| remains almost insensitive to the changes in G5o (Figs. 2a, 
2b and 3c). 

The results summarized above are robust. Indeed, the same tendencies in horizontal fluxes 
are observed in 2D vs. ID horizontal variabilities, in presence vs. absence of vertical 
inhomogeneity, in variable vs. flat cloud top and base, with strong surface albedo vs. a “black” 
surface, and, finally, with Henyey-Greenstein vs. realistic phase function with the same 
asymmetry parameter. Of course, the magnitude of horizontal fluxes increases with the gradients 
in cloud structure. Stronger cloud variability, clear sky gaps in cloud structure and oblique 
illumination may substantially enhance the magnitude of horizontal fluxes as additional sources for 
photon horizontal transport. All these factors make the situation very complex and often when 
analyzing high resolution satellite images one cannot distinguish between the competing effects of 
radiative smoothing and geometrical shadowing for both transparent and absorbing channels. 
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Figure Captions 

Figure 1: Dependence of pixel-by-pixel IP A accuracy II///HI on averaging scale r and single - 
scattering albedo G5 0 . Horizontal distribution of cloud optical depth is simulated with bounded 
cascade models (p = 0.3, H = 1/3, <T> = 13, 10 cascades, pixel size = 25 m). Flat cloud top and 
cloud base, geometrical thickness h = 300 m. No cloud gaps are added. Henyey-Greenstein phase 
functions with asymmetry parameter g = 0.85 is used. Surface is “black.” The results are 
averaged over ten independent realizations. Solar zenith angle 0o=O°, single-scattering albedo GJ 0 

- 1.0, 0.999, 0.99, 0.98, 0.95, 0.9. (a) Reflectance, ll/fflll; (b) Transmittance, II// 7 -II; (c) 
Absorptance, II// 4 II. 

Figure 2: Pixel-by-pixel IP A accuracy for reflectance, transmittance, absorptance, their sum, 
and the total horizontal fluxes H. Cloud model and scattering conditions are the same as in Fig. 1 . 
(a) e 0 =0°; (b) 0 o =6O° 

Figure 3: Wavenumber spectra E(k) for the IP A and MC calculated radiation fields. Cloud 
model and scattering conditions are the same as in Fig. 1 . Three single-scattering albedos G 3 0 are 
used: 0.999, 0.98, and 0.95. (a) 0 o =O°, nadir radiance; (b) 0 o =6O°, nadir radiance; (c) 0 o =O°, 
transmittance (cloud optical depth is added for reference). 













Fig. 2a 



Fig. 2b 
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